Trans-Planckian dispersion and scale-invariance of inflationary perturbations 
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We investigate the insensitivity of the predictions of inflationary models with respect to mod- 
ifications of Planck energy physics. The modification we consider consists in replacing the usual 
dispersion relation by nonlinear ones. This way of addressing the problem has recently recieved 
attention and contradictory results were found. Our main result is to show that the adiabaticity 
of the mode propagation and the separation of two scales of interest, the Planck scale and the 
cosmological horizon scale, are sufficient conditions for the predictions to be unchanged. We then 
show that almost all models satisfy the first condition if the second is met. Therefore the introduc- 
tion of a nonlinear dispersion is unlikely to have any discernable effects on the power spectrum of 
cosmological perturbations. 
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A number of recent publications |l], g, [|, || || address 
the so-called trans-Planckian problem of inflationary cos- 
mology: in some inflationary scenarios, scales of cosmo- 
logical interest today were redshifted from wavelengths 
far below the Planck length, / p i , at the onset of in- 
flation. Since we do not have at present a consistent 
quantum theory at these scales, this appeal to unknown 
physics questions the validity of the standard predicitions 
of inflation, i.e. scale invariance of the power spectrum of 
cosmic fluctuations and Gaussianity of their distribution. 

In usual scenarios, the metric fluctuations responsible 
for the generation of large scale structure and the temper- 
ature anisotropies of the cosmic microwave background 
radiation (CMBR) originate from quantum fluctuations 
of a massless field propagating in a background geome- 
try whose expansion rate H is determined by the poten- 
tial energy of the inflaton field H [7|. More specifically, 
the power spectrum of the metric fluctuations can be 
obtained from the vacuum expectation value of the two- 
point function of a minimaly coupled scalar field, eval- 
uated at the time when the field wavelength is equal to 
the cosmological horizon H" 1 [||. If the potential is suf- 
ficiently flat, H is nearly constant in time, giving rise 
to a nearly scale invariant horizon crossing amplitude. 
Moreover, the linearity of the field equation results in 
the Gaussianity of the distribution statistics when the 
state of the field is vacuum. 

To question the robustness of these results with respect 
to changes of Planck energy physics, Refs. [|[ [| [| mod- 
ified the propagation of the quantum field in this high 
frequency domain, in close analogy to what was done for 
Hawking radiation (see Q for a review and references). 
The modification consists in replacing the linear disper- 
sion relation which follows from Lorentz invariance by ad 
hoc nonlinear dispersion relations. The linearity of the 
field equation was maintained, guaranteeing that Gaus- 



sianity be preserved. Different forms of the dispersion re- 
lation were investigated, all of which reducing to the lin- 
ear behavior for small frequencies but becoming nonlin- 
ear above a given wavenumber n c which can be thought 
to be of the order of the inverse Planck length. Assuming 
that horizon crossing occurs in the linear regime of the 
dispersion relation (i.e., that H <C k c ), Refs. || || find 
no violation of scale invariance of the fluctuation power 
spectrum. On the other hand, Ref. j|] claims to find 
significant scale non-invariance in some cases. 

These opposite results raise three questions we shall 
address in this Letter: (1) What are the sufficient re- 
quirements for a nonlinear dispersion relation not to sig- 
nificantly violate the standard predictions of inflation? 
(2) How do the first deviations from the standard results 
manifest themselves in the power spectrum ? (3) Is the 
cosmological situation different from the black hole case 
H [h], [Tl|, [l2| where no deviation in the thermal spectrum 
was found when the cutoff frequency is much higher than 
the Hawking temperature, and if yes, why? 

Given the spatial homogeneity of the background ge- 
ometry, the propagation of a linear field is governed by 
the time evolution of the field modes Xk with conserved 
comoving wavenumber k. Therefore the fc-th component 
of the power spectrum is determined only by the ampli- 
tude of Xk at the (conformal) time r] k of horizon crossing, 
i.e., when k — a(r/k) H, where a(r/) denotes the scale fac- 
tor of the Friedmann-Robertson- Walker metric. When 
using the conventional dispersion relation, Xk obeys the 
equation 



Xk + u 2 Xk = o 



with u> 2 given by 
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where dashed variables stand for derivatives with respect 
to 77 (hence, H = a' /a 2 ). Upon modifying the dispersion 
relation, the wave equation is still given by Eq. ([!]) but 



to 2 is replaced by 



uj 2 F (v) = Hv)F(k/a( V ))f 



a"(v) 
0(17) 



(3) 



F(k/a) is an as yet unspecified function which converges 
to the usual physical wavenumber k = k/a for k <C k c . 

The phase and amplitude of \k are fixed by two condi- 
tions. First, its normalisation is provided by the Wron- 
skian condition, 



XkXk* ~ Xk Xk 



(4) 



which follows from the requirement that the creation and 
annihilation operators satisfy the canonical commutation 
relations [ jl3| . This condition equally applies for uiq and 
luf- Secondly, one needs a positive frequency condi- 
tion. Then the mode so defined uniquely specifies the 
creation operator at that can be used to construct the 
/c-component of the field state in Fock space. 

The physical requirement fixing the state at rji, i.e., 
at the onset of inflation, is that the fc-components corre- 
sponding to physical wavenumbers K\ = k/a(rji) that are 
very large, i.e., well above the temperature and H , should 
be in their ground state. When using the usual frequency 
luo of Eq. (0), there is no problem in charcaterizing the 
vacuum since ujq is constant for early 77;. However, ujf 
might still be time dependent at rji. Hence one looses 
the conventional way of defining the vacuum. In order 
to define it, we shall follow the conventional adiabatic 
approach E3], In an adiabatic regime, the adiabatic vac- 
uum equivalcntly follows from the mathematical proper- 
ties of the modes Xk [ H or from the physical require- 
ments to minimize either the probability for excitation 
of a comoving particle detector or the energy density of 
the field || . This is because both quantities are based on 
the local properties of Xk that are accurately described 
by WKB waves. 

It is well known that the adiabatic vacuum is associ- 
ated with the positive frequency WKB solution of Eq. ([!]) , 
normalized according to Eq. (0) and given by 
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The key point for us is that the probability for non- 
adiabatic transitions fbif (i.e., transitions leading to spon- 
taneous excitation of the vacuum and interpreted in sec- 
ond quantization as pair creation) is determined by the 
corrections to Eq. (||). Indeed, in the adiabatic regime, 
i.e. when the adiabatic parameter 



(6) 



satisfies C(ry) <C 1 for all 77, < 77 < 77/ , the pair cre- 
ation probability amplitude at 77/ is equal to backscat- 
tering amplitude (3 k , i.e. the amplitude of the negative 
frequency part of Xk evaluated at 77/, when starting at 
r\i with the positive frequency mode Eq. (||) . Explicitely, 
Pk is defined by decomposing the exact solution \k at 77/ 
as 



Xk(Vf) = a kX r*(rif)+f3k{xr B (Vf)} 



(7) 



where \a k \ 2 - \(3 k \ 2 = 1. 

In the context of inflation, if the evolution from 77/ to 
the horizon crossing time 77^ is governed by ujq (see below) 
[3 k modifies the power spectrum (A$)^ ~ k 3 \xk(Vk)\ 2 /a 2 ■ 
Using Eq. (Q), we find 

\Xk(Vk)\ 2 = \xf\vk)\ 2 [1 + 29?e(/? fc e*-) + 0(f3 2 )} (8) 

where Xk ^ s ^ ne un perturbed amplitude evaluated using 
luo when working in the unperturbed vacuum ((3k = 0). 
e l< ^ k is a pure phase which depends on the value of 77/. 
Hence (3k determines the amplitude of the modulation 
around the usual horizon crossing amplitude. 

The condition for Eq. (^) to be valid is that the horizon 
and cutoff scales be well separated: 



mpi 
H 



» 1 



(9) 



In this case, 77/ can be chosen such that the physical 
wavenumber k = k/a is, at the same time, much smaller 
than k c and much bigger than H. In this interme- 
diate regime, the second term in Eqs. @ and (||) is 
small compared to the first for t] < r/f and can there- 
fore be neglected in the evaluation of C. Furthermore, 
F is assumed to be in the linear regime at 77/, thus 
F(k/a(rj > 77/)) = k/a. The evolution of Xk from 77/ to 
77fc is therefore identical to the standard, non-dispersive 
case. Consequently, to obtain identical horizon crossing 
amplitude, it suffices to show that the modified mode 
governed by lof coincides at 77/ with the purely positive 
frequency mode obtained when using loq. This requires 
only (3 k < 1. 

In the case where, in addition to C <§; 1, C — > for both 
77 — > r\i and 77 — * 77/, (3k is exponentially small fTij . In our 
specific case governed by ujf this exponential would be 
of the type e~ Kc ^ H , i.e. completely negligible if Eq. (||) 
is satisfied. However, the condition C — > for 77 — > 77/ is 
unlikely to be met since one has to evaluate (3k at a given 
finite time, in particular before the last term of Eq. (||) 
becomes significant. When evaluating (3k with C(rjj) ^ 0, 
only weaker predictions can be made. Nevertheless, when 
C € 1 and slowly varying from 77^ till 77/, (3k is bounded 
by the maximum of C(r/). (This can be obtained from 
the boundary contributions at 77/ generated by multiple 
integrations by parts, see Sect. 2.4 in |l4[|.) Hence (3k is 
still negligible when C(r]) <C 1. 
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In brief, the robustness of the preditions of inflationary 
models follows from (a) the linear relationship between 
the modulations of the horizon crossing amplitude and 
the backscattering amplitude and (b) from the smallness 
of this amplitude when both adiabaticity, C « 1, and 
scale separation, a 3> 1, are met. 

This is in complete agreement with what has been 
found in black hole physics. In that case, due to the sta- 
tionarity of the metric, one works with fixed frequency to 
rather than fixed wavenumber k and the differential anal- 
ysis proceeds with the radial variable r rather than with 
the time ij. When the appropriate adiabatic condition is 
met, see Eq. (47) in ficfi , and when the cutoff scale is well 
separated from the Hawking temperature, which is fixed 
by the inverse radius of the hole and which plays the role 
of H, positive norm WKB waves provide good approx- 
imations of the modified modes, in the sense that the 
amplitude of negative norm component of these modes 
is negligible. This guarantees that the modifications of 
the Hawking flux will be equally negligible. In fact, the 
first order modulations of the flux have been shown to 
scale like er~ 2 Thus, as in cosmology, scale sepa- 

ration and adiabaticity of the mode propagation protect 
the predictions from alterations generated by a nonlinear 
dispersion relation. 

Let us now show that in most cases, scale separation 
Eq. (0) leads to adiabaticity through the relation C ~ 



<C 1. To this end, we write Eq. (||) as 



C 



H 
~F 



HndF 
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(10) 



with n — k/a. Obviously, the dispersion relation used 
by Unruh |l which behaves as F(n) — > k c for k — > oo 
is globally adiabatic as the first term is ~ er -1 «C 1 and 
the second term vanishes. This is in agreement with the 
results of Refs. § |. 

We consider next the class of functions F(k) which 
grow monotonically with k. Here, the first term in 
Eq. ([l0|) can be neglected as it is smaller than the cor- 
responding term obtained using the linear dispersion re- 
lation and which is already negligible. The second term 
also vanishes asymptotically for F ~ K n with any positive 
n, as well as for exponentially and logarithmically grow- 
ing functions. This class includes the dispersion relation 
of Ref. H in agreement with the author's conclusion. On 
the other hand, it also includes the superluminal disper- 
sion relation investigated in Ref. Q where it is claimed 
that it gives rise to deviations from the standard result. 
Our result is obviously in conflict with that of Ref. . 

Finally, the dispersion relation may become singular 
at a finite wavenumber 



F(k) = 



(1 ~ k/Koo)" 



(11) 



For n = 1 , the second term in Eq. (|l0|) becomes equal to 
the first at k — n,^, whereas it vanishes for n > 1 and 



diverges for n < 1. Only in the latter case, the notion 
of an adiabatic vacuum becomes invalid near the cutoff 
scale, making the predictions for inflationary fluctuations 
dependent on an alternative choice of initial conditions, 
or a later specification of the initial state, see below. 

The situation is different if F decreases monotonically 
for k > Ko ~ k c , as the first term in Eq. ( |To| ) grows with- 
out bound and C ^ for nonlinear F. Even for strictly 
positive F, i.e. disregarding the additional complications 
due to complex dispersion relations, there exists a time r\ 
where F is smaller than H. According to Eq. (0), this 
causes the breakdown of the adiabatic condition because 
cosmological expansion begins to dominate the dynamics. 
In other words, those modes with wavenumber 3> H^ 1 
are "frozen" in the beginning 1 . As they are streched by 
cosmic expansion, their proper frequency grows accord- 
ing to the nonlinear behavior of the dispersion relation. 
Therefore, because of Eq. (JoJ) , the adiabaticity condition 
will be satisfied before the dispersion relation reaches the 
linear branch. As an example, consider the subluminal 
dispersion relation of Corley & Jacobson [|ll] (also ana- 
lyzed in Ref. f§): 



F(k) 



M 1- 



(12) 



At k = kq/2, the corresponding adiabaticity parameter 
is C ~ (J- 1 < 1 for n = 0(1). 

If an early non-adiabatic phase is followed by an inter- 
mediate adiabatic regime on Planckian but sub-horizon 
scales, any viable initial solution cannot deviate too far 
from the adiabatic vacuum without violating the consis- 
tency of the inflationary framework. In the intermediate 
adiabatic regime, the number of these primordial parti- 
cles is constant and the usual rules for computing the 
energy density apply. Therefore their energy density will 
be of the order of n^nik) where n(k) — \(3k\ 2 is their oc- 
cupation number density. Without some fine tuning of 
n(k), this term may become comparable to or even larger 
than the inflaton potential energy, invalidating the key 
condition for inflation to take place. This point was re- 
cently (and independently) developed by Tanaka fl6|| 2 . 
On the other hand, if the contribution of n(k) to the 
expansion rate is negligible compared to that of the in- 
flaton potential, one has n(k) <C c~ 2 -C 1 for k/a ~ k c , 
i.e., these modes must be extremely close to the adia- 
batic vacuum. Thus the predictions of inflation remain 



1 This phenomenon was also investigated in Ref. Jli| which ap- 
peared shortly before this article was finished. 

2 For an early discussion of this point in the black hole context 
we refer to hit. In this case, for exactly the same reasons, the 
hypothesis of having a regular metric near the horizon severly 
constrains the occupation number n(u) if an intermediate adia- 
batic regime exists. 



4 



unchanged. These remarks apply, for instance, to the 
case n < 1 in Eq. ( pi] ) if k^/H <C 1 and to subluminal 
dispersion relations with F < H . 

To summarize, we found that all regular, monotoni- 
cally growing dispersion relations are globally adiabatic 
and thus unlikely to give rise to non-standard perturba- 
tion spectra if horizon crossing occurs in the linear dis- 
persion regime. Only those functions that decrease with 
increasing wavenumber, or that contain a certain type of 
singularity, violate adiabaticity at early times. In those 
cases, the specification of the initial state is conceptu- 
ally less clear than in the former where the adiabatic 
vacuum is available. However, if adiabaticity is restored 
well before horizon crossing, the initial conditions can be 
assigned in this intermediate regime, guaranteeing that 
the adiabatic solution holds after this point. Again, the 
standard inflationary predictions remain unchanged. 

One can also question the validity of the assumption 
of scale separation itself, and speculate about potential 
signatures if the field is still in the (slightly) nonlinear 
dipersion regime at horizon crossing. Given that in the 
simplest scenarios for inflation, the CMBR fluctuation 
amplitude scales inversely with a evaluated at horizon 
crossing of the respective wavelength , we know that 
a w 10 5 . Whether or not this number should be consid- 
ered "large" cannot be answered without independent ex- 
perimental constraints on residual Planck scale effects on 
long wavelength phenomena. Indeed, Amelino-Camelia 
& Piran |l8| invoke nonlinear corrections to the disper- 
sion relation as a possible explanation for the observa- 
tion of ultra-high energy cosmic rays and TeV-photons 
whose fluxes would normally be suppressed by interac- 
tions with CMBR and infrared background photons, re- 
spectively. Intriguingly, their allowed region of parameter 
space is consistent with a second order correction whose 
coefficient corresponds to a ss 10 5 . Nonlinear disper- 
sion acting on length scales comparable to the horizon 
length during inflation therefore cannot be ruled out on 
experimental grounds alone. By noting that the ampli- 
tude of the WKB solution Eq. ([|) depends only on the 
local value of F rather than its history, one can estimate 
the ratio of the horizon crossing amplitude in the stan- 
dard and dispersive theories as ~ H/F(H) || (assuming 
strict adiabaticity of the prior evolution). Like in the 
standard theory, time translation symmetry of de Sit- 
ter space has to be broken by making H a weakly time 
dependent variable in order to obtain a scale dependent 
perturbation spectrum. Depending on g — AF(k)/Ak, 
with Ak = Ak a covering the cosmological scales observ- 
able today, we speculate that the natural reddening of the 



spectrum due to the slow decline of H will be slightly en- 
hanced (g > 1), reduced (0 < g < 1), or even reversed 
(g < 0) as the fluctuation amplitude changes more or 
less rapidly, or in the opposite direction, as a function of 
proper wavelength than in the standard theory. 

In conclusion, nonlinear dispersion alone is unlikely to 
have any clearly discernable effects on the power spec- 
trum of cosmological perturbations. As suggested in 
Ref. H and argued on the basis of noncommmutative ge- 
ometry in Ref. |l9| ] , a more promising signature of Planck 
scale physics may be non-Gaussianity of the fluctuation 
statistics. 
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